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Abstract 

We determine the fermionic corrections to the nucleation rate of bubbles at the 
electroweak phase transition. The fermion determinant is evaluated exactly and 
by using the gradient expansion. The gradient expansion is found to be a reliable 
approximation and is used to extrapolate to the large values of v n = (2n + l)nT 
needed in the Matsubara sum. The contribution to effective action is found to be 
negative and to be given, essentially, by the gradient terms, the finite part of the 
wave function renormalization. Only the top quark contribution is evaluated, it 
is of the same order as the Higgs- and W-boson contributions found previously, 
but of opposite sign. 
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1 Introduction 



The physics of the electroweak phase transition has been discussed recently in 
various aspects [[I]]. Many subjects, as e.g. the question of baryogenesis |2|, |3| are 
still controversial ||, |5|. Even the nature of the phase transition is not known at 
present. The temperature dependence of the effective potential has been studied 
in perturbation theory |], [7| as well as in lattice simulations |8|, ||. If the mass 
of the Higgs is not too high (less than My/) the phase transition is supposed 



to be first order |]I0[ . In this case the transition from the symmetric vacuum 
with massless particles to the broken symmetry phase would proceed via bubble 
nucleation. This phenomenon as well as its cosmological aspects, has been studied 



by various groups jO], 0' 0- 



Part of the basic information needed in developing the bubble formation and 
expansion scenario is the determination of their nucleation rate. In the small 
temperature span of 1 GeV in which the phase transition takes place bubbles 
of various sizes can be formed. Their nucleation rate varies over several orders 
of magnitude. The basic rate is determined [14j 15, |16| by the classical minimal 
bubble action. Its exponential is the tunneling rate. The semiclassical reaction 
rate includes, however, also preexponential factors, the fluctuation determinants, 
determined by the fluctuation of W boson, Higgs and fermion fields in the back- 
ground of the minimal bubble profile. The bosonic fluctuations have been com- 
puted recently [|T7|, [TB| and found to yield sizeable suppression factors, the one- 
loop effective action (or equivalently the free energy divided by the temperature) 
being of the same order as the classical bubble action. In the high temperature 
theory, obtained by retaining only the Matsubara frequency 0, fermions do not 
contribute. However, recent determinations of the fermionic contribution to the 



sphaleron rate [19[] let us expect that at least the top quark will influence the 



transition rate in an essential way. Of course we have to perform such a com- 
putation in the four- dimensional finite temperature quantum field theory, i.e. by 
summing over all Matsubara frequencies. 

The plan of this paper is as follows: in the next section we will introduce the 
model and set up the basic relations for the bubble nucleation rate. The fermion 
determinant is defined in section 3. In section 4 we will discuss the leading terms, 
of first and second order in the external field vertex function. Part of their contri- 
bution is contained already in the effective potential and should not be included 
again. This renders the discussion in that section very lengthy and pedantic. In 
section 5 we present the computation of the finite higher order contributions 
to fermionic effective action. The results of an exact numerical computation are 
compared with an analytic approximation based on the gradient expansion. The 
latter one is used, then, to obtain the actual results which are presented and 
discussed in section 6. 
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2 Basic relations 



As a basis for the computation of fluctuations we need an action which describes 
the nucleation of minimal bubbles and which can be used to compute the bubble 
profiles. This must necessarily be a finite temperature action which describes 
the essential features of the phase transition. We will be using here an action 
obtained in the electroweak theory by evaluating the one-loop effective potential 
a la Coleman- Weinberg |2(| . This has been computed by various authors f2l| , p2j 
23]; here we use the formulation of Dine et al. The finite temperature action 



was found by these authors - for the temperatures relevant to the phase transition 
- to be well represented by its high temperature limit. Though we will compute 
the fermion determinant at finite temperature, it should be a good approximation 
to compute the bubble profile from this approximate action the mores so as 
the one-loop corrections will be computed exactly and their contribution to the 
approximate one-loop effective potential will be subtracted at the end in order to 
avoid double-counting (see below). This way of determining the 'classical' profile 
from an action including already one-loop corrections has been discussed critically 
by Weinberg f2M; we think that this intermediate approach nevertheless yields 



relevant information, namely the order of magnitude and sign of the corrections 
as well as their structure and relative importance. 

The basic Euclidean finite temperature action then given by | 



P 



o 



Sr. = I dr I d 3 x 



+ S F . (2.1) 



$ is the complex doublet of Higgs fields. Here this field will always occur as a 
background field describing the minimal bubble. It can be parametrized then as 

<£(£) = v $(x) (J J . (2.2) 

Vm is the high temperature potential which includes the one-loop effective po- 
tential in the high temperature approximation 

VLt(* + *) = D(T 2 - T 2 )&$ - £T$ f $ 3/2 + ^$t$2 _ 

Its parameters are given - for 6 m = - by 
D = (3m 2 w + 2m 2 )/8v 2 
E = 3g 3 /32n 
B = 3(3rr4 - 4m 4 )/64vr 2 t; 4 

T 2 = (m 2 H -8v 2 B)/AD (2.4) 

2 2 
771 771 

X T = A-3(3m 4 w ln-^-4m^ln-^)/167r\ 4 (2.5) 



3 We leave out the gauge fields here, since they do not appear in the bubble background 
field configuration and since we are not going to compute their fluctuation determinant. 
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with In = 21n47r — 27 and In ap = 2hi7r — 27. 

For T > T the potential has a minimum at |<I>| = corresponding to the 
symmetric phase and a second minimum at 



l»l=«XT> = ^ + V(^)W) P.6) 

where 

, 2 (t) = ^(T 2 -T 2 ). (2.7) 

This minimum is degenerate with the one at $ = at a temperature defined 
implicitly by 



T c = To/y/1 - E*/D\ To . (2.8) 
Tc marks the onset of bubble formation by thermal barrier transition. 

The fermion action Sf can be written - for vanishing gauge fields and for 
Higgs field configurations of the form fl2.2|) - in terms of four-component Dirac 
spinors as 



Sf = f dr I d 3 x 



. / //' 



(2.9) 



where the Yukawa couplings are related to the fermion mass matrix via g Y Vq = 
vrv* . The sum over / is over flavors and colors. Here we will consider only the 
contribution of the top quark. This has been done already in the high temperature 
action given above for the reason that its contribution is much larger than the 
one of lighter quark and lepton fields. This will also be the case for the exact 
one-loop action. 



The process of bubble nucleation is - within the approach of Langer |M and 
Coleman and Callan [15, 16], followed by the work of Affleck [^5|, Linde [26| and 
others - described by the rate 



Y/V 




exp(— S) 



Jf\ 1/2 
Jb) 



(2.10) 



Here S is the high-temperature action, Eq. ( |2.1| ), with the new rescaling, min- 
imized by a classical r independent minimal bubble configuration (see below). 
Jfib are the fermionic and bosonic fluctuation determinants which describe the 
next-to-leading part of the semiclassical approach. Jp - whose computation is the 
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aim of this work - will be defined below; its logarithm is related to the fermionic 
one-loop effective action by 



°eff 



hxJ F 



(2.11) 



Finally U- is the absolute value of the unstable mode frequency. 

The classical bubble configuration is described by a vanishing gauge field and a 
real r independent spherically symmetric Higgs field. For the bubble configuration 
we make the Ansatz 

= u(T)0(r) ( J ) (2.12) 

where we have rescaled Eq. ( |2.2j ) via t>o$(r) = v(T)(j)(r). In all our numerical 
computations we will use the scale (g^t^T)) -1 for the coordinates. Defining the 
rescaled high temperature potential as 



la 



V 



2)<? 



with 



AET 



X T v(T) 

and the high temperature coupling 



4 
3 



1 - 



v{T)\ 



the bubble action is given by 

§ 4w 



r dr 



v(T) 



1 d0 

2 [fa 



(2.13) 



(2.14) 



(2.15) 



+ V ht (<f>) 



(2.16) 



~gl{T) Jo 

It is minimized if <f>(r) is a solution of the associated Euler-Lagrange equation 

- 0"(r) - V(r) + = (2.17) 
r d(p{r) 

with the boundary conditions 

lim 0(r) = and 0'(O) = . (2.18) 

r— »oo 

The bubble configuration varies from small thick wall bubbles to large thin wall 
bubbles in the narrow range (~ 1 GeV) between T and Tc, both of order 100 
GeV. We will use the variable 



y = 3(l-e/2), < y < 1 
instead of T to parametrize this range of temperatures. 



(2.19) 
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3 The fermionic fluctuation determinant 

The fermionic action Sf of quarks can be rewritten in four component Dirac 
notation and for time-independent background configurations as 



S F = J o drj ^( T ^-m(x))^ (3.1) 

13 dr f d 3 xip\d/dr - H)tp (3.2) 
u J 

where m(x) = gyv(T)(j)(r) and 

H = 7 Q (-r/V + m(x)) . (3.3) 

The field fluctuations are subject to antiperiodic boundary conditions at r = 
and t = j3, i.e. 

ijj(x,(3) = -^(x,0) (3.4) 

which determines their frequencies to be the Matsubara frequencies v n = (2n + 
1)ttT with integer n, — oo < n < oo. Integrating out the fermion field leads to 
the fermionic prefactor 

.1/2 _ Una(Wn + Ug 

where u> a denotes the eigenvalues of H. Using the fact that these eigenvalues 
occur in pairs ±cu a we can rewrite this as 



,7-1/2 llnaK^n -r iu Q , . 

<?F = r, 777. , ,QN I 3 " 5 ) 



The fermion contribution to the effective action is therefore given by 

^>*Hl lndet (^)- (37) 

Here u n = (2n + l)irT, the fluctuation operators M. and Ai° are defined as 

M = H 2 = -A + V(x) 

M° = (H°) 2 = -A (3.8) 

and 

V={ -^f^f) ) . (3.9) 

[ iaVm{x) m (x) J 

A method for computing such fluctuation determinants numerically has been 
described recently [ 27| . Before we discuss the numerical part of the computation 
we have to ensure that the quantities we are going to compute are finite. The 
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effective action as defined formally in Eq. ( |3.7| ) is divergent. This is easily seen by 
expanding it w. r. t. the potential V. One generates then the series of Feynman 
graphs depicted in Fig. 1, of which the first and second order graph are divergent. 
We will see that our numerical method allows to separate these two graphs from 
the remaining series which can be computed exactly and is finite. The divergences 
of the two leading graphs are obviously those of ordinary perturbation theory. 
Their divergent parts can be cancelled by the counterterms of the T = theory. 



This will be decribed in the next section, following the work of Refs. |28fl and 



4 Renormalization 

We note first that the effective potential is renormalized usually at T = in 
such a way that the vacuum expectation value vq and the Higgs mass, i.e. the 
second derivative at the broken symmetry minimum, are kept at their physical 
values. This will determine the counterterms in the effective action. We will, on 
the other hand, expand the effective action around the symmetric vacuum defined 
by $ = 0. The discussion of renormalization - especially at finite temperature - is 
therefore somewhat cumbersome; we will follow all steps very explicitly, though 
this may seem somewhat pedantic. 

In order to cover both kinds of expansion we rewrite the effective action as 

S8e _--l„ l7 _--j : ^ detM _ A + (4) J (4.1) 

where now fip — rrip = gyVo if one expands around the T = broken symmetry 
vacuum and /xp = if one expands around the high temperature symmetric 
phase. The 4x4 matrix U is then given by 



u = v-A = \ m ^\~-f 1'"f V * !■ ■ < 4 2 > 



The field $ is normalized in such a way that it takes the value 1 if the Higgs field 
is at its T = vacuum expectation value. 

In the following we will discuss the contributions of first and second order in 
the external potential both at T = with massive propagators and at large T 
with massless propagators. Anticipating that we work at finite temperature we 
replace the space-time integration by (3 J d 3 x. Furthermore, for the second order 
diagrams we separate space and loop momentum integration formally, though 
the external momentum q acts as an operator in x space. This compact notation 
should not lead to confusion. 
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For the first order diagram at T = we have 

SH\*,0) = -\W J d 3 x{m 2 F \<S>{x)\ 2 J ■ (4-3) 

Using a 4 - momentum cutoff A we get 

S$(*,0) = ~P J d 3 x(m||$(x)| 2 -^)(A 2 -^ln(A 2 /^)) . (4.4) 

Before we dispose further by introducing appropriate counterterms we turn to 
the second order contribution at T = 0. We have to evaluate 

Sg^O) = Uf3 J d 3 x{{m 2 F \^{x)\ 2 -^ 2 F f + m 2 F \V^{x)\ 2 ) (4.5) 
r d 4 p 1 1 

x J (2tt)V + /4 ( p + q y + ^ F ■ 

A standard computation, using again a 4-momentum cutoff A leads to 

S$(M = ^Jd 3 x((m 2 F mx)\ 2 -^ F ) 2 + m 2 F \W^\ 2 ) (4.6) 



x 



(In - 1 - f 1 da ln(l + ^a(l - a))) 

\n F J Jo n F 



Putting this together with the first order contribution and requiring all corrections 
to the potential and the wave function renormalization to vanish at q 2 = we 
find that we need a counterterm 

A 2 



c.t. 

eff 



8tt 2 



(3 J d 3 x(m 2 F (\<5>(x)\ 2 -l)) (4.7) 
^(ln - l)P J d 3 x(mUmx)\ 4 - 1) + m 2 F \V® (x)\ 2 ) 



^m F (3 J d 3 x{\$l 



x)\ 2 ~l) 



8tt 2 

If the first and second order diagram are evaluated in the symmetric vacuum we 
find 

A 2 r o ,,o 1 /A 2 



x(3 J d 3 x(m F \<Z>(x)\ 4 + m F \V$(x)\ 2 ) . 

Leaving out the field-independent terms - which are due to a difference in the 
vacuum energy density - we get 

S^ +2 ' S \$,0) + S C J- = -^(ln^-2)/3/rf 3 xKl^)| 4 + ^|V$(x)| 2 ) 



which is to be evaluated with the background profile. The last term is already 
contained in the high temperature potential ( [2.3|) , it occurs in the coefficient 
—DTq = —(rrijj — 8vqB)/A via the top contribution to B. 

Turning now to the first order diagram at finite temperature we have to eval- 
uate 



d 3 p 



1 n i-oo „ 

(4.10) 

In order to separate T = and finite temperature contributions we use |28|| 

OO 

T „£x,((2n+l)*rT) 2 + P 2 + /£ (4 ' n) 

1 1 1 



2E F E F exp(E F /T) + 1 

with Ep = Jp 2 + ftp. In the last line the second term vanishes as T — > 0, 
so the first one represents the T = contribution which we have considered 
earlier. Inserting the second part into the expression for J 1 ^ we find the finite 
temperature part 

AS«(*,T) = 2/J / d 3 x(™ F Mx)\ 2 - fj.%) / t^-t ■ (4.12) 

J J (2ti) 6 E f exp(E F /T) + 1 

We have to evaluate this expression for a bubble in the symmetric vacuum where 
$ 7^ only locally and where fi F = at spatial infinity. Then E F = |p| and we 
find 

AS£($,T) = ^m 2 F (3 J d 3 xmx)\ 2 . (4.13) 

This contribution is already taken into account in the T 2 term of the three di- 
mensional high temperature action (p.l[) . Therefore the finite temperature part 
of the first order tadpole diagram has to be omitted entirely. 

As the second order contribution at finite T we have to evaluate 

S$($,T) = -P J d 3 x4(m F \<5>(x)\ 4 + m F \V$(x)\ 2 ) (4.14) 



+oo 
n=— oo 



d 3 p 



(27r) 3 (^ + p2)( v 2 + (p + q) 2 ) 



Momentum integration and Matsubara frequency summation can be carried out 
via 



T E 



d 3 p 1 



(27r)S(^ + p 2 )(i^ + (p + q) 



S 



+00 „J . 

T ^L da I 



d 3 p 



n=— oo 



(2vr) 3 (z/2 + P 2 + a (i - a ) q 2)2 



+5? ri rj^^d 1 

^ '» 7 (2vr) 3 d{p 2 ) vl + p 2 + ail - a)q 2 " 1 J 



n=— oo 



1 /-i , /■ d 3 p -d 7T 2 T 

da - — — -. n . — — tanh(-£/Q,/2T) 
in y 2vr 3 d n 2 2£„ V 7 7 



7r 2 Ty 7 (2vr) 3 d(p 2 ) 2E Q 

where i? 2 = p 2 + a(l — a)q 2 . The T = part may be recovered by performing the 
limit T — > 0. Subtracting this part we find for the finite temperature supplement 

1 lo^ dpp2 2ptp^ a eMEjT) + l = (416) 

1 , -1 f 00 dp 1 
da- 



An 2 Jo E a exp(E a /T) + 1 



so that 



AS%($,T) = -^/3|d 3 x(m||$(x)| 4 + m||V$(x)| 2 ) (4.17) 



fl /-oo 

x / aa 



lo E a exp(E a /T) + 1) ' 

Part of this term is already contained in the high temperature effective action. 
The momentum integral has been considered by Dolan and Jackiw [^8|; it can be 
expanded to leading order in T, i.e. up to terms of order 1/T 2 as 

If this is inserted into the previous equation we find 

l(2) ^ _ 1 f 1 j fa(l-a)q 2 



A^(*,T) = Y^J o da(\n ^ j + Try) (4.19) 

(3 J d 3 x{m F \<5>{x)\ 4 + m 2 F \V${x)\ 2 ) . 
The integration over a can then be performed, replacing the first parenthesis by 

( ln (^) " 2 + 27) • (4.20) 

If this is added to the zero temperature result Q4.9|) the term \n(q 2 /m 2 F ) — 2 in 
this equation gets replaced by ln(T 2 /m 2 F ) + 2 ln7r — 27. This term appears in the 
high temperature potential in At- 
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Collecting from Eqs. ( f4.9| ) and ( |4.17| ) the terms which have not yet been 
included into the high temperature potential we find the following renormalized 
contribution of the first and second order Feynman graphs: 

A^ n + 2 )($,T) = f^J\^)\ 2 + ^Mv)\ 2 ) (4-21) 

f 1 , f°° dp 
da 




E a exp(E a /T) 



Here we have taken into account the color factor 3, we have replaced mp by m ( 
and we have replaced the formal x integration by a q integration in order to deal 
with the q dependence of the kernels correctly. Note that $ 4 (x) becomes here the 
square of the Fourier transform of $ 2 . 



5 Calculation of the finite part of the effective 
action 

What remains to be evaluated is the sum of all Feynman graphs of order 3 and 
higher. This contribution is finite, we denote it as S^(4>). We present an exact 



numerical computation, using a general theorem on functional determinants |29 
and an analytic approximation based on the gradient expansion. The results of 
the two methods will be compared at the end of this section. 

5.1 Numerical computation 

As mentioned above the fermionic one-loop effective action at finite temperature 
can be written - including the color factor 3 - as 

^ eff = -^ £ mdet =-3X>^M- (5.1) 

z n=-oo \ n J n=0 

As the background field is spherically symmetric the determinant can be decom- 
posed into its partial wave contributions. This is readily done by introducing the 
usual spinors for given j and I = j ± 1/2 as given in the textbook of Bjorken and 



Drell [301. One finds then 



s^-EE 6 ( 2/ + 1) ^Jiiyn) ■ (5.2) 

n=0 Z=0 



Here the the partial wave determinants Jiiv) are defined as 



JiW) = det 



' Mi + v 



2 
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with the partial wave fluctuation operators 



Mj = Mj' + Vtr) 

Mf = -xf^UlH': 1 ) °A ,5.3) 



dr 2 rdr I r 2 1 (/ + !)(/ + 2 



I m 2 (r) dm(r)/dr 1 , N 

V r = { j / \ / j 2 / \ \ • (5.4) 

1 dm(r)/dr m (r) 1 

A very fast method for computing fluctuation determinants is based on a theorem 



on functional determinants [|29[ which can be generalized to a coupled (n x n) 
system: 

Let f(V, r) and f°(z/, r) denote the (n x n) matrices formed by n linearly 
independent solutions ff{y, r) and /f°(z/, r) of 

(M I + z/ 2 ) y j;(z/,r) = (5.5) 

and 

(M° + z, 2 ) y jfW) = (5.6) 

respectively, with regular boundary conditions at r = 0. Here the latin lower 
index denotes the n components while the different solutions are labelled by the 
Greek upper index. Let these solutions be normalized in such a way that 

limf(z/,r) (f°(z/,r)) _1 = 1. (5.7) 

r— >0 V / 

Then the statement of the theorem is: 

det (Mi + z/ 2 ) _ detf(z/,r) 
Jl[V) = det (MP + v 2 ) ~ ™ detf°(i/,r)- (5 - 8) 

Here the determinants on the left-hand side are determinants in functional space 
while those on the right-hand side are ordinary determinants of the nxn matrices 
defined above. The theorem already has been applied to the calculation of the 



one-loop effective action of a single scalar field on a bubble background in |p7 
and of a fermion system at temperature T = on a similar background in |31 
which we refer to for more technical details. 



In the numerical application the solutions f£ are written as |32 



fk(v,r) = (% + ht(v,r))i li (vr) (5.9) 

with the boundary condition h% (u, r) — > as r —>■ 0. Of course the value 4 depends 
on the channel. This way one generates a set of linearly independent solutions 
which near r = behave like the free solution as required by the theorem which 
then takes the form 

J{u)= lirndet(^ + /^,r)). (5.10) 
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The functions h%(v,r) satisfy the differential equation 



ar z \r % i k \y r ) ) dr 



v,r) = V w {r)(6% + h%( V ,r) 



ilAvr) ' 



(5.11) 

This equation can easily be used for generating the functions hf order by order 
in the potential V. Introducing the contribution of order k in the potential as 
and and defining via 

oo 

h ( k) _ j2 

j=k 



as m 



31], the relevant contribution S^ s is found to be 



where 



^ = -E^ (3) K 



n=0 



(5.12) 



/C®( 



LXJ 

5^6(2i + l) Hm {indet (l + h (1) (zv 



1=0 



tr ( h (1) (i/, r) + h (2) (i/, r) - - [h (1) (v, rf 2 



(5.13) 



The expression K,^\v) was evaluated by summing the partial waves computed 
numerically up to l max = 30 extrapolated for higher values of I using an Ansatz 
al~ 5 + bl~ 6 + cl~ 7 . The asymptotic behaviour is supposed to set in at values of 
/ >> vR where R is the typical radius of the bubble. Since R has typical values of 
20-40, this means that for our l max the extrapolation becomes unreliable already 
for v of order 1. We will discuss this point again below. 



5.2 Analytic approximation using the gradient expansion 

In JJTJ an approximation of the gradient expansion type has been given for the 
one-loop effective action at zero temperature for the case of a massive fermion 
with Yukawa coupling to an external scalar field. The calculation to be done here 
is similar because the Matsubara frequency u n enters in the same way as a mass 
term and so the structure of the determinant is of the same type. 
The logarithm of the determinant J(v) can be written exactly as 

\nj(v) = Tr ln(l + G V) (5.14) 

where the free Green function Go is defined by 

(M + v 2 ) G = 1. (5.15) 
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Defining the Fourier transform of the potential V as 

V(q) = J d 3 x V(x) exp(— iqx.) 
the fluctuation determinant can be expanded as 



(5.16) 



(ftp Jy f d 3 qj 



V(q,-) 



\nj(u) =g tr — ^ y (27r)3 iy ^ (p + Q . )2 + 



:(2vr) 3 5( 3 )(Q A 



(5.17) 



where 



Qi = Eq*- 



Expanding the denominators including terms up to order v 2k 4 we get 
k 1 1 

2 (p + Qi )2 + ^2 - ( p 2 + V 2) 



k q2 

h (p 2 + 



+ 



3 (p 2 + V 



2\2 



EQA+EQ' 



1=1 



(5.1* 



As the potential is V = m F $ + 771^7 V$ we have 



tr HV(qj) ~ 4m^ 

1 



(5.19) 



After inserting these expansions in Q5.17T ) and transforming back to x-space the 
remaining integrations (except of one space integration) and the fc-summation 
can be done. Of course we have to omit those terms which are divergent. As we 
are working in four dimensions these are the terms with k = 1 and k = 2, they 
have been discussed in the previous section. 

Our final result written in the form of /C(z/) = ?>Yn.J'{y) is 



m^(V$) 2 
2z7 



4 
3 



+ 





m F § 2 


V 2 




m 2 F $ 


V 2 



1 - - 



3m^$ 2 3m^$ 4 



2 v 2 



m F § 2 



(5- 



3z/ 2 



1-1 



m 2 F <fr T 



This expression has to be evaluated using for $ the numerical bubble profiles. 
The accuracy is only limited by the accuracy of these profiles and by that of the 
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numerical integration. With our numerical precision the results are reliable to at 
least 6 significant digits. 

The approximate results for A3 3 ) can be compared with the exact numeri- 
cal ones computed using (|5.13|) . For the purpose of comparison we treat v as a 
continuous parameter. We display the exact and approximate results in Figs. 1 
and 2 for two typical bubble profiles, a small bubble with y = 0.6 and a large 
bubble corresponding to y = 0.3 (see ( |2.19| ) for the definition of y). The analytic 
approximation is seen to describe the trend of the exact results over the whole 
range. The gradient expansion is expected to converge at large v. This expecta- 
tion is substantiated by the exact numerical results in the region where they are 
reliable. It is seen, however, that the exact results start dropping off at values of 
v ~ 1 — 2; as mentioned above this is related to the fact that the convergence of 
the partial wave summation becomes poorer with increasing v. Since the values 
of v relevant for the Matsubara frequency summation ( |5.12| ) are v > ttT « 9 ( in 
our units gv) we have to rely on the gradient expansion in computing the finite 
temperature effective action. 



6 Results 

We have computed the finite temperature fermionic effective action for Higgs 
masses of 60, 70 and 80 GeV and for top quark masses m t = 160, 170 and 180 
GeV. These results are given in Tables 1 to 3. As mentioned in section 2 we 
have considered only the contribution of the top quark since lighter quarks and 
leptons will give negligible contributions. Their contribution has already been 
dismissed in the basic high temperature action ( [2.1|) and including them would 
be inconsistent. 

For each set of mass parameters we have determined the bubble profiles for 
various values of the variable y defined in (|2.19| ); y determines the temperature 



and the bubble action. We give separately the renormalized first and second 
order contributions AS^ 2 \ Eq. (|4.21|) determined in section 3 and the finite 

sum of all higher order contributions S^ s whose computation was discussed in 
the previous section. It is given by the analytic expression Eq. (|5.20| ), inserted 
into the Matsubara sum (|5.12|) . The total one- loop effective action, reduced by 



the terms included already in the 'classical' effective potential is given, of course, 

by the sum AS^ +2) + S$ . 

The fermion determinant is seen to yield a negative contribution to the effec- 
tive action, which means that bubble nucleation is enhanced by this contribution 
(cf. Eqs. ( ^.1(J| ) and ( 2.11| )). It is interesting to analyze the relative importance of 



the various contributions. The (3) part given in Eq. ( |5.20| ) is relatively small, also 



the zero-derivative part of the leading order (1 + 2) contributions in Eq. (|4.21 



This is not unexpected, as the zero- derivative part was already included into the 
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effective potential, so its leading first and second order contributions are sub- 
tracted. The dominant contribution is the gradient term in Eq. (|4.21|) , the finite 
part of the wave function renormalization. It is also displayed in Tables 1 to 3 in 
the colums labelled A^f\ This means that the top quark contribution can be 
described, essentially, by local terms in analytic form: those already contained in 
the effective potential and the finite part of the wave function renormalization; 
this is very convenient because these terms can be incorporated into the basic 
action from which the bubble profile is computed. If this is done, the remaining 
corrections can be expected to be very small. In a selfconsistent determination 
of the bubble profiles such an action can be expected to yield an excellent first 
approximation. 
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Table Captions 



Table 1 Corrections to the effective action at finite temperature for m# = 60 
GeV. e and y are defined in Eq. ( gig ) and ( gig ); A^+ 2 ) is_ defined by Eqs. 
( |4.21|) , ASj^ffl is the gradient part as explained in section 6; S'® is given by Eqs. 



(|5.12[) and (|5.20| ); S, the classical action of Eq. ( [2.161 ) is given for comparison. 



Table 2 Corrections to the effective action at finite temperature for m# = 70 
GeV. Definitions as in Table 1. 

Table 3 Corrections to the effective action at finite temperature for m# = 80 
GeV. Definitions as in Table 1. 



Figure Captions 

Fig. 1 The loop expansion of the effective action. The lines represent the prop- 
agators and the dots indicate the vertices V(x). 

Fig. 2 Results of the numerical computation compared to the analytic approx- 
imation for y = 0.3, mjj = 60 GeV and m t = 170 GeV. The dots interpolated 
by a dotted line represent the numerical results, the solid line is the analytic 
approximation of Eq. fl5.13| ) 

Fig. 3 The same as Figure 2 for y = 0.6. 
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m t [GeV] 


T [GeV] 


e 


y 


AiS (l+2) 


A o(l+2) 


o(3) 
D cfT 


S 


160 


94.894 


1.867 


0.2 


-23.160 


-23.250 


-18.923 


308.02 




94.854 


1.800 


0.3 


-9.637 


-9.753 


-3.774 


132.30 




94.682 


1.600 


0.6 


-1.685 


-1.693 


-0.0781 


24.847 


170 


94.5288 


1.867 


0.2 


-23.752 


-23.695 


-20.418 


278.47 




94.495 


1.800 


0.3 


-9.973 


-10.091 


-4.184 


121.40 




94.347 


1.600 


0.6 


-1.734 


-1.742 


-0.0854 


22.686 


180 


94.6605 


1.867 


0.2 


-24.558 


-24.291 


-20.929 


251.60 




94.631 


1.800 


0.3 


-10.002 


- 10.114 


-4.110 


107.27 




94.5056 


1.600 


0.6 


-1.757 


-1.766 


-0.0855 


20.245 



Table 1 



m t [GeV] 


T [GeV] 


e 


y 


AiS (l+2) 
^ren 


A c(!+ 2 ) 
*^grad 


o(3) 
D cff 


S 


160 


106.015 


1.867 


0.2 


-18.933 


- 18.740 


-6.746 


217.14 




105.980 


1.800 


0.3 


-7.956 


- 8.006 


-1.384 


94.690 




105.826 


1.600 


0.6 


-1.387 


-1.391 


-0.0282 


17.630 


170 


104.7850 


1.867 


0.2 


-20.091 


-19.899 


-8.412 


207.39 




104.7545 


1.800 


0.3 


-8.438 


-8.500 


-1.726 


90.411 




104.6202 


1.600 


0.6 


-1.470 


-1.475 


-0.0352 


16.830 


180 


104.094 


1.867 


0.2 


-20.897 


-20.743 


-9.602 


193.72 




104.0674 


1.800 


0.3 


-8.751 


-8.799 


-1.946 


83.921 




103.951 


1.600 


0.6 


-1.522 


-1.528 


-0.0396 


15.632 



Table 2 



m t [GeV] 


T [GeV] 


e 


y 


A £(l+2) 
rcn 


A o(i+2) 


o(3) 
°cff 


S 


160 


117.541 


1.867 


0.2 


-15.433 


-15.177 


-2.512 


156.45 




117.510 


1.800 


0.3 


-6.483 


-6.483 


-0.514 


68.050 




117.374 


1.600 


0.6 


-1.136 


-1.139 


-0.0107 


12.708 


170 


115.4848 


1.867 


0.2 


-16.747 


-16.411 


-3.388 


153.28 




115.4573 


1.800 


0.3 


-6.985 


-6.984 


-0.688 


66.460 




115.337 


1.600 


0.6 


-1.222 


-1.225 


-0.0142 


12.410 


180 


114.0042 


1.867 


0.2 


-17.696 


-17.287 


-4.163 


146.18 




113.9801 


1.800 


0.3 


-7.387 


-7.386 


-0.852 


63.649 




113.874 


1.600 


0.6 


-1.287 


-1.290 


-0.0174 


11.846 



Table 3 
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